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I. INTRODUCTION 



The modification of fundamental laws of physics with the change of observational length 
scale is the subject of the renormalization group (RG) [1]. Through the RG flow equation one 
may probe the dependence of the effective coupling constants on the characteristic length. The 
otherwise complicated flow pattern becomes rather simple in the vicinity of the fixed points where 
the linearized RG flow along with scaling provide a recipe for classifying the coupling constants 
via their dependence on the characteristic scale. The irrelevant coupling constants are those 
which decrease as the scale is moved toward the infrared (IR) direction. The physical content 
of the theory is insensitive to the actual choice of these coupling constants in the IR end of 
the region where the linearization of the RG equation is applicable. Within this regime where 
the usual concept of universality is recovered, the physics is parameterized by the others only, 
namely, the relevant and marginal coupling constants. 

In the realistic models we find several scaling regimes when the renormalized trajectory 
passes by different fixed points. There are fixed points for the Theory of Everything, GUT, 
Standard model, QCD and QED, to mention some of them. Although the true renormalized 
trajectory approaches all of them for certain values of the cutoff, it reaches the first one only. 
In fact, in the scaling regime of, say the fixed point of QCD, some of the interactions of the 
Standard Model generate nonrenormalizable vertices in terms of the quark and gluon fields [2]. 
These vertices deflect the renormalized trajectory from the fixed point as we move up in energy. 
(Hence it is physically not too crucial whether or not the ultraviolet (UV) fixed point really 
exists. All we shall assume here is the scaling up to a certain energy scale.) The traditional 
goal of local field theory is to give an account of these vertices in terms of elementary particle 
exchanges in a manner which is renormalizable at higher energy. Yet renormalization, i.e., the 
removal of the UV cutoff is necessary only for the Theory of Everything. In fact, when the scaling 
is investigated at the other fixed points then the higher energy reactions always make these fixed 
points unstable in the UV direction. 

The scenario sketched above leads to a serious limitation of the use of the concept of uni- 
versality. It is true that there are "islands" of autonomous scaling regimes where physics can 
be parameterized by relevant or marginal operators only, but these operators usually vary with 
scaling regimes and the matching of relevant coupling constants is rather nontrivial. Even though 
we can establish the the importance of certain coupling constants in a given energy range, the 
physics at a different scale will be governed by different set of coupling constants [3] . 

There is one last fixed point as we move towards longer distance scale, the infrared (IR) 
fixed point. Macroscopic physics is characterized by the scaling at this IR fixed point. Can this 
fixed point have relevant operators ? The answer is negative for theories with a mass gap. To see 
this it is sufficient to recall that the dependence of the coupling constants on the cutoff is to take 
into account the effects of the modes which are eliminated as the cutoff energy is lowered. At 
energy well below the mass gap the fluctuations are suppressed and the evolution of the coupling 
constants slows down. Thus the IR limit of the theory is stable. 

The situation is more interesting for theories without mass gap. Realistic theories with 
spontaneous symmetry breaking belong to this class. For such theories, IR divergences can pile 
up and generate relevant coupling constant in the IR regime. It is the main result of this paper 
that this indeed happens in simple four dimensional scalar models. In this case the long distance 
physics of the model is not universal, i.e., it cannot be parameterized completely by the relevant 
and marginal coupling constants of the UV fixed point. 
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In arriving at this result one needs an improved version of the RG equations. The usual 
method of renormalizing the theory is to follow the mixing and the evolution of selectively few 
coupling constants. One traditionally chooses those coupling constants which are relevant or 
marginal in the vicinity of the UV fixed point. As we lower the cutoff to the natural mass scale 
m,R of the theory, the scaling properties change fundamentally. In the IR side of the natural 
mass scale where k < uir, one deduces the scaling laws corresponding to the IR fixed point. 
However, the UV and the IR scaling regimes are separated by a crossover at ttir, and there is no 
reason whatsoever to expect the same set of relevant or marginal coupling operators for both UV 
and IR fixed points. Furthermore we do not have a simple power counting argument to find out 
the relevant operators of the IR fixed point. Thus in order to establish the scaling operators of 
the IR fixed point, we have to trace down the evolution and the mixing of many more operators 
which might well be irrelevant in the UV scaling regime. This can be achieved with a RG flow 
equation which is capable of handling the mixing between infinitely many coupling constants. 
Such an improved RG equation has been obtained in [4] and was subsequently applied for the 
UV scaling regime in [5] in the leading-order approximation of the derivative expansion for the 
renormalized action. We present in this paper the RG equation [4] applied in the next order of 
the derivative expansion. This allows us to verify our claim about the existence of the relevant 
operators at the IR fixed point in the first two orders of the derivative expansion. 

The organization of the paper is the following: In Sec. II we give a brief derivation of the RG 
equation [4] in the leading order of the derivative expansion for in the one-component scalar 

model and show the emergence of the IR singularities in certain (3 functions. Section III contains 
the technical details for deriving for the O(N) Xcp 4 theory with two distinct wave function 

renormalization constants, Z^^ and Z^j for the longitudinal and the transverse components, 
respectively. In Sec. IV, we derive a set of three coupled nonlinear RG flow equations for 
Zk,e. and Zk jt . Asymptotic scalings in both UV and IR limits are discussed in Sec. V. Section 
VI contains our conclusions. Two appendices are included supplementing the detail of deriving 
the RG equations for the paper. 



II. ONE-COMPONENT SCALAR FIELD THEORY 



Our starting point for deriving the RG equation for a system characterized by the field (f>(x) 
is to introduce the coarse-grained "block variable" : 



v 



<Pk(x)= I p k (x-y)(f>(y), (2.1) 

where 



J = J d d x = fl, (2.2) 



via a smearing function pk(x), with k 1 being the characteristic linear dimension of the region 
over which the field averaging is performed. In this paper, we shall choose pk(x) to be 

(2 - 3) 
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or pkip) = Q(k —p), i.e., a sharp cutoff [6]. Although pk(x) acts as an upper cutoff, we shall 
use A as the /c-independent UV cutoff for the theory. 

Given a set of blocked variables 4>k{x), the blocked action Sk can be deduced from 

e -s fc (*) = I D[(j)] Y[5((/> k (x) -$(x))e- s ^\ (2.4) 

J X 

where the field average $ of a given block is chosen to coincide with the slowly varying back- 
ground. By performing the functional integration in loop expansion subject to the 5 function 
constraints one finds 



Q(<b\-Q(<b) Itt'i &2s I [' dS f 8 2 S \ OS 

k{ } ~ U + 2 %(p)00(-p) 2 J p d<l>{p) {d<P(p)d<P(-p) ) dcP(-p) 

= S($) + -Tr'lnK - - / FK~ X F, 
2 2 J 



(2.5) 



where F and K are the first and the second functional derivative of the bare lagrangian, respec- 
tively, and Tr' denotes the trace sum over internal space as well as the restricted momentum space 
with k < p < A. How blocking transformation modifies the propagator A (a; — y) = K~ x (x — y) 
can be seen explicitly by considering a free scalar theory: 

f F(p)A(p)F(-p) — f F(p)A(p)F(-p), (2.7) 

where 

A(x-y)= / ——e(p-k) (2.8) 

Jp P Tf 

is the "blocked" propagator with an effective IR cutoff scale k. In the limit k — > 0, one recovers 
the original A(x — y). 

Eq. (2.5) is far too complicated so the derivative expansion [7] is used at this point. The 

form 

W) = f; / d d xL%\<S>{x)) (2.9) 

n=0 J 

is assumed where L^\<&(x)) is a homogeneous polynomial of order 2n in the space-time deriva- 
tives. We shall truncate the expansion at n = 1, retaining only the wave function renormalization 
function Z k {$) and the blocked potential £/fc(<E>). Such truncation, being justifiable in the IR 
limit for high enough space-time dimension d, yields simpler differential equations when substi- 
tuted into (2.5). In principle, however, equations that generate the scale dependence of Zk(Q) 
and the higher order derivative terms must also be calculated in the framework of the derivative 
expansion in order to have a closed system. 
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It is worth mentioning that (2.5) gives the one-loop effective potential for k = [6]. The 
modes with nonvanishing wave number are eliminated independently in the one-loop approxima- 
tion. However, the result can be greatly improved by the successive elimination of the degrees 
of freedom. In the improved scheme, the contribution of a particular mode which has been 
integrated out is kept for the elimination of the next mode, thereby taking into account the 
interactions between the modes. By decreasing the cutoff infinitesimally from k — > k — Ak. we 
generate from (2.5) the evolution of the potential, = £4(3?), which for $(x) = $ and 

d = 4, becomes: 



U k _ Ak ($) = U k ($) + Ak^\n 



Z k (<$>)k 2 + d%U k {$) 
Z k (0)k 2 + dlu k (o) J 



(2.10) 



or equivalently, a differential equation of the form: 



kd k U k ($) 



k 4 



16tt 2 



rln 



Z k ($)k 2 + d 2 U k ($) i 
Z k (0)k 2 + d*U k (0) \ 



(2.11) 



in the limit Ak — > 0. This equation describes the renormalization of the potential with arbitrary 
dependence on the field $. The solution for U k =o(&) differs from the usual one-loop effective 
potential mentioned before insofar that the effects of the operators which are irrelevant at the 
UV fixed point are retained during the elimination of the degrees of freedom. This difference is 
negligible for a weakly coupled theory as long as no new relevant operators are generated outside 
the UV scaling regime. 

Equation (2.11) can be derived by resumming the one- loop contributions to evolution of the 
potential. In fact, it can be written as 



kd k U k {$>) = 



k 4 



16tt 2 

k 4 
16tt 2 



In 



In 



-Z k ($)k 2 + m 2 k 




Z k (0)k 2 + m\ J 



+ In 



+ 



l+iZkWtf + ml^dlVkW 



dlV k {<$>) 



Z k ($)k 2 + m 2 



+ 



where d%V k {$) = d%U k {§) - ml and m\ 



d%U k (0). 



Z k ($)k 2 + m 2 

(2.12) 

The last line contains the sum of the 
one-loop graphs with increasing number of d\V k {&) insertions. The external legs of these graphs 
which are attached to $ in 9|,Vfc($) are carrying zero momentum and the modes with momentum 
k are propagating along the loop. This is just the set of graphs one has to sum up in eliminating 
the modes with momentum k. 

There are certainly higher loop corrections to (2.12). However, the terms of the order m in 
loops contains m integrations over a <i-dimensional shell in the momentum space. When only few 
modes are eliminated, Ak ~ 0, the integration over each shell yields, on the dimensional ground, 
a new small parameter: 

C = x> ^ 

which helps suppress the higher loop contributions to the RG evolution equation in (2.11). This 
is the basis of the "exactness" of for a RG equation which is formally obtained in the one-loop 
approximation. 
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As the IR limit is approached with k — > 0, ( becomes ill-defined. To examine the behavior 
of ( in this regime, we introduce an IR cutoff by considering the system in a box of size L. With 
the number of degrees of freedom being N d = {^) d , the momentum integral measure takes on 
the form: 




(2.14) 



in term of the dimensionless momentum, p^. In the "ultimate RG transformation", where only a 
single mode is eliminated each time, the small parameter would be -j^, the inverse of the number 
of degrees of freedom remained. Therefore, the IR limit can be reached in two different manners: 
(i) The limit L — > oo is taken first before k — > 0; and (ii) k — > ^ is first taken for finite systems 
followed by L — > oo. We immediately notice that £ can be kept small only for case (i) and 
becomes 0(1) for case (ii). In another words, the limits k — > and L — > oo do not necessarily 
commute. In fact, the gap of the two dimensional cr-model, whose existence is established with 
reasonable accuracy using procedure (ii) is absent if (i) is employed instead [8]. 

It seems that if the two limits k — > and L — > oo are not commuting, procedure (ii) would 
be more reasonable to describe the dynamics of local interactions. In that case the RG equation, 
(2.11), can only represent a partial resummation of the perturbation expansion. Such a loss of 
the effectivity is due to the presence of the length scales, L and /c -1 , in a system without IR 
mass gap. Since in the subsequent treatments we implement the one-loop RG equation for the 
IR regime, the conclusions drawn from such computation are strictly relevant only for case (i). 
It remains to be seen if they can be carried over to the procedure (ii). 

Adopting (i) as our approach, we find indication of the emergence of relevant operators at 
the IR fixed point for massless theories due to the following handwaving argument: For the sake 
of simplicity we keep the wave function renormalization constant, Zfc($), to be unity in (2.11) 
and introduce the ft functions for the coupling constants for $ n as: 

ft n (k) = d%kd k U k ($), (2.15) 

which can be obtained by substituting (2.11) into (2.15). These ft functions describe the evolution 
of the coupling strengths of small fluctuations around the constant background <&(x) = Let 
us now consider a model where 

<9tt/ fc =o(0) = 0. (2.16) 

This is what one commonly calls a "massless" theory since U k =o($>) is just the effective potential. 
However, this name is misleading when spontaneous symmetry breaking occurs with < <p{x) >^ 
since d%Uk=o(< 4>{ x ) >) is now non-vanishing. The leading IR contribution for the (2.15) comes 
from the highest power of k 2 + <9|[/fc(0) in the denominator: 

which shows that for d$U k (0) > the coupling constants for the odd powers of the field blow 
up in the IR limit. Although we have not found the scaling operators there ought to be relevant 
ones which drive the IR divergences. 
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Such conclusion could have been reached by considering the graphs which contribute to the 
evolution equation: 



16n 2 *\Z k ($)k 2 + m 2 2\Z k ($)k 2 + m\ 



In the UV scaling regime where k 2 » 9|L r fc($), the dominant contribution comes from the 
graphs with least number of propagators. The evolution of the vertex with n legs is described by 
joining two legs of the (n + 2)th order vertex having momentum k. One reproduces the usual (3 
functions, e.g,: 

for the mass squared, and in the usual one-loop approximation to the <^ 4 model where there is 
no sixth-order vertex at the cutoff, 

v ; 16tt 2 (k 2 + m|) 2 

in the next-to-leading-order approximation in (2.18). 

As we enter the IR regime with k 2 « <9| the scaling laws quickly change. For a theory 
with mass gap the contributions are O(^r) and the evolution slows down indicating the absence 
of relevant operators. But for massless theories, lim/ c _ > o TO l = 0, the dominant contributions 
are received from graphs with the maximal number of propagators between the vertices. In the 
absence of other dimensional parameter, we find <9§Vfc(0) ~ k 2 and the IR contribution to the 
evolution of the n-th order vertex is dominated by the one-loop graph with n insertion of the 
vertex <9§ as shown in (2.17). 

Unfortunately this result is not interesting. The theory develops a nonvanishing vacuum 
expectation value for the field either due to the masslessness or the presence of the odd powers of 
the field in the potential. Thus, the coupling constants computed at $ = are not characterizing 
the strength of the interactions of small fluctuations in the vacuum. The true vacuum with 
< <&(x) >y^ shields the IR divergences. However, if the theory possesses a continuous symmetry 
which is broken spontaneously then Goldstone's theorem guarantees the presence of the massless 
modes in the vacuum. The more careful repetition of this simple argument for the iV-component 
scalar field theory is the subject of this work. 

III. DERIVATIVE EXPANSION FOR THE O(N) MODEL 

We consider a generalized bare O(N) scalar field lagrangian of the form 

C(<P) = ±Z a (d^ a ) 2 + V(<j ) ), (3.1) 

with a = 1, • • •, N. The theory is chosen to be in the symmetry broken phase and the the 
direction a = 1 is chosen to be in the expectation value of the field. The extra subscript in Z is 
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used differentiate between two wavefunction renormalization constants, one for the longitudinal 
component and the other for the transverse ones: 

2^, a = 1 

Z a = { (3.2) 
Z u a = 2,---,N. 

Alternatively, one may write (3.1) as 

£(0) = ^Md^ 1 ) 2 + \z t {d^f + V(cf>) 
= - 1 -Z a <f> a d 2 <f> a + V{<j>) 



(3.3) 



via the relation: 

d 



dcj) 1 

Z t = -^(Z t cf> i ), i = 2,---,N. 



(3.4) 



We shall split <j>(x) into the slowly varying background x( x ), an d the fast-fluctuating modes £(x) 
such that 

r x a (p), o< P <k 
<p a ( P ) = { (3.5) 

{ £ a (p), k < p < A. 

Noting that 4>k(p) = Pk{p)4>{p) = x{p)i one then integrates out the fast-fluctuating modes £(x) 
by using the loop expansion to obtain the blocked action Sk(&) as a function of the blocked field 
average $ = ($ 2 ) 1 / 2 : 

5 fc (*) = -ln f D[ X ]D[Z}l[5(Mx) ~ *(*)) exp{-S( X + 0} 

= -ln /" DMjJafoCp)-^)) 
^ p 

x | exp{-S( X ) - fc(p)F a (-p) - \Je{p)K ah {p, -p)?{-p) + ■ • } 
= -ln | D[ x ]]|]5( x (p)-$(p))exp{-5(x)-^ r IV / lnK+i^ FK~ 1 F + • • • j 

= + ^Tr' In K($) - 1 y 

= 5($)+^($)+^ 2 ($), 

(3.6) 

where 

F«($) = ^ = -i^ a) ($)$ c a 2 $ c - ^Z a ($)(<9 2 $ a + $ a d 2 ) +F (a) ($), (3.7) 



K ab {§) 



d 2 s 

l ) a (x)d(j) b (y) 



= j-iz( a&) ($)$ c <9 2 $ c - Z a {<$>)5 ab d 2 + V {ab \$) 

- -z[ a \<s>)(d 2 <s> b + $ 6 a 2 ) - -z( b X$)(d 2 § a + $ a a 2 ) \5 4 (x - y ), 

2 2 j 



(3.8) 



and 



an y 

v ' d<f) ai d<f) a2 ■ ■ - d(j) a ™ 

y(aia 2 —a n ) _ - ; 



(3.9) 



d(j) a ^d(j) a2 ■ ■■d(j) a " 

Note that Tr' denotes the trace sum over the internal symmetry space as well as the restricted 
space-time. 

As noted in the Introduction the blocked action can be expanded as 



(3.10) 



For the computation of Zk(§) and Uk(&) it is best to choose a non-constant, slowly varying 
blocked field which is written as 

$ a ( x ) = <$>% + 4> a ( x ), (3.11) 
with $g = f&o^' 1 - Simple comparison of (3.10) and (3.6) gives 



(3.12) 
(3.13) 



and 

£ a ,fc($o) = Z a (Q ). 

From now on, quantities with no written arguments are understood to be evaluated at $o- 

In order to obtain higher order correction for the wavefunction renormalization constant, we 
incorporate the effect of 4> up to quadratic order by writing 

R ab = ( R ab + dR ab + dR ab + ^6^4^ _ y ) + q^3 ^ ^3) 

F a = F a + 6F? + 5F* + O(0 3 , d 3 ), 



(3.14) 



where 



o 

a 6 



K% b = -(Z a 5 ab + Z^^qS^S^O 2 + V {ab) 
SKZ b = ^o(2^1 a) 5 a ' 1 5 b ' 1 - Z ( b a) S b ^ - Z^S a '^d 2 

-Z[ c ^ c 5 ab d 2 --(z { b a) d 2 4> b + z^d 2 r) - \zf h) ^d 2 ~^ 

2 V / 2 

-- \z { b ac) <p c (d 2 4> b + fid 2 ) + z^frid 2 ^ + ^ a d 2 ) + z^frd 



5Kf = 



(3.15) 



5Kf 



2lc 



+ Z { a cd ^ c ^ d 5 ab d 2 ] + -V( abcd ^ c d , 

2 
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and 



F a = -Iz a $ <5 a,1 <9 2 + y (a) 



1 



SF? = ~ Z { b a) ^ S bA d 2 (f) b + Z a (8 2 a + 4> a d 2 ) + Z^^o^ftd 2 + V (a 



z ( b a) 4> b d 2 4> b + z( a V<f> s c ' 1 j) b d 2 4> c + z^4> b (d 2 4> a + 4> a d 2 ) 



(3.16) 



1 



+ Z { a bc) ^ 5 aA ^ b ^ c d 2 + -V {abc) (f) b (f) c . 

Note that summation over repeated indices (a,b,c,d= 1, • ■ • , N; i, j,k,= 2, • • • , N) is implied, 
and the numerical subscripts correspond to the order of (p. The effective action can now be 
written as: 

S fc ($) = S($ + 4>) + 6Sl($o +4>) + 6S 2 ($ + $) (3.17) 

where 

5($ o + 0)= J (-^-0 a d 2 (P a + V + V( a) a + ^4> a V^4> b + •••), (3.18) 
SSl = ^Tr'ln^o + SK + 5K X + 5K 2 ^ = ^Tr'ln^o + SKq) - |Tr' (k^SKxKq-HK^) 

1 r / \ n -i 1 r / \ n 

+ -^(-l)-Tr'[(K - 1 ^ ) K _1 ^iJ + 2E(- 1 ) nTV lro" 1 ^o) ^ _1 ^2 



n=0 



n=0 



+ 



(3.19) 



and 



SS 2 k 



= \ I (F + SFi + 5F 2 ) (K Q + 5K + 5K X + 5K 2 ) _1 (F + 5F X + 5F 2 ) + 
2 J p 

= i jf jFoKo-^l + Ko-^^-^o + ^o^^l + ^o" 1 ^)" 1 ^ 

+ ^Kq-^I + K^SK^Fq - FqK^SK^I + K^5K )- 2 F 

+ F K^(1 + K- 1 5K )- 1 5F 2 + 6F 2 K^(1 + K^SKq^Fq 

- FoK^dK^l + K^ 1 SK )- 1 5F 1 - SF 1 K^ 2 SK 1 (1 + K^SKq^Fq 



(3.20) 



+ F Ku L [(Kq '8X^(1 + Kq L 5K q )- 6 - (Ko L K 2 )(l + K^ l 5K ) 
+ SF 1 Kq 1 (1 + K^ 1 5K )- 1 5F 1 ^ + ■■■. 

To illustrate the above formalism, we consider the O(N) scalar Xcp 4 theory defined by 

A 



V(cf>) = ^ 2 {x) + ^ 2 {x)) 2 , 



(3.21) 
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where 
With 



2 = 



V™ = (p 2 + ^ 2 )5 ab + ^S^S 1 ' 1 , 

o 



(3.22) 
(3.23) 



the matrix K$ h in the momentum space representation has the eigenvalues 

Aj 1 = Z^p 2 + lie, a = 1 
A^ = Z t p 2 + u u a = 2,---,N, 



(3.24) 



U£ 



^ 2 + 



(3.25) 



which yields 



Note that (3.4) implies 



Tr'lni^ 6 = Q J I InAj 1 + (N - l)lnA t _1 I. 



(3.26) 



z, 


= Z e + Z i £ 1) $ 




= 2ZP + Z^ 


< zr 


= 3Z^ + ■■■ 




= z f + z£ i) $ 




= z? j) + ..., 




z t = z t 



(3.27) 



(i) 



(3.28) 



Z + W = 2Z, (l) . 



In evaluating (3.19) and (3.20), we act on the re-dependent <j> in the trace by the derivative 
operators contained in F and K. The trace is performed in the plane-wave basis by commuting 
the momentum operator p^ = id^ to the right end of the expressions. The commutation relations 
utilized for this procedure are tabulated in Appendix A. Note that whenever an operator p M acts 
on (p it yields id^cf) which contributes in the IR because its Fourier decomposition is vanishing 
above the scale k. When the operator p^ reaches the right end of the expression it gets re- 
placed by the trace integration variable, p M . Although this momentum value is in the ultraviolet, 
the contribution represents a simple number which multiplies the <p dependence in the infrared. 
Therefore, a matrix element of the form TV \F((p)p^Q((p) can be separated into Tr' ' J-{4>)id^{Q{4>)) 
and Tr / jF(0)^((/))p M using the commutation techniques. Iterating this algorithm, the contribu- 
tions to Z which are of the form Tt'J r ((p)d 2 Q((p) can be isolated. Moreover, since the blocked 
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action is real we can actually commute the derivative operators to the left end instead within 
the trace, as was done in [7] and adopted here in this paper. 
One arrives after lengthy algebra at 



Z^S ab + ^o^V' 1 + Z^S a ^)U c + \z^$> + \z^ 



Z <f)8 ah + -^ (Z { b ac) 5 b ^ + Z^S^ild 2 ^ - \zi a) d 2 ft + Z^d 2 ^ a 
2 . 

- \z™* Q &p - 2i[z^5 ab + \Mzi ac) ^ 1 + Z^5^)]p,d^ c 

+ V ^4> c - i Plt [zl a) d^ b + zf>a/] + • • • , 
- - [z { b ac) 4> c d 2 4> b + z^4> c d 2 r + z^frd 2 ^] + ■■■, 

2 

6F? = -p 2 \z a (f) a + Z^QoftS"' 1 ] + 5 ab u t ^ b + ^QffiS"' 1 - Z a d 2 (f) a 



_ I (zl a) 5 h ^ + Zf^"' 1 ) ^d 2 ~4> b - ip„ [z a d^ a + zl b) Q 6 a >% 



+ 



and 



(3.29) 



(3.30) 



(3.31) 



5F? = l -p 2 [Z< 6 >#° + zf c) ^8 a ^ b - \ [Zj)j« + Z^p + Zf b) ^]d 2 ^' 

+ lv {abc) 4> b 4> c + ---. 

2 



(3.32) 



As for the calculation of 8S%, since no x integration occurs in (3.20), it is only the Fourier 
transform of </>. If (p(p) is constrained such that p < k, then the contribution of <5S| to the 
blocked action vanishes since the p integration is performed over the range k > p. 

In order to identify the contribution to the wavefunction renormalization constants and the 
blocked potential we now rewrite the blocked action £&(<&) as: 



\ztrd 2 ^ b + u k + (u { k a \r) + u ( k ab) p) + 



(3.33) 



from which one obtains the (^-independent blocked potential: 

U k = V+^Tv'\n^K + 5Koj 

l 

= V + U P ( lnA 7 1 + ( Ar - 1 ) lnA *" 1 )- 
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(3.34) 



For the wavefunction renormalization constants, we first collect terms proportional to ^d 2 ^ 1 for 
computing Zk,i in the longitudinal direction. Use of (3.27) and (3.28) yields (see Appendix B): 



Zk,t = Ze + ^ / jaii + a'u + an + a* ± - 26n| 



1 



A 2 



Z e u e A}\ 2 <f> 2 - Z £ (1) A$ (1 - u e At + 2u\A\) - (Z £ (1) ) p 2 (l + u\A\) 



+ (N - l)u t A 3 t { A t [ZtY^o ~ Vu t ] + ^o4 1} (1 " 2«tA t ) } 
- zf ^ A?(ziV + ^o) - i(Z,W) VA,A t (4 + rz?A? + n 2 A 2 ) 

+ 3Z / (11) #A, + zf ) A t 



(3.35) 

One may replace the Z factors above by the corresponding Z's since the difference involves 
contributions having more than two orders of derivative in the Z's. Finally, to obtain Zk,t(&), we 
simply replace $o in Zk,i(®o) by <E>, with implied "normal ordering" such that all p dependences 
be moved to the front of the ^-dependent expressions. As noted in [7], there is no ambiguity 
in this procedure provided that we carefully compare the terms in the expansions of (3.10) and 
(3.17). 

In a similar manner, we can write down the ^-dependent wavefunction renormalization 
constant in the transverse ij direction as: 



Z t + 



a*j -26ij} 



= Z t + - 



-Zjf p 



u 



2 Aj + A 2 (3 + 



2N + 3 



-u 



2 A?) 



(zf >) WA?(2 + u 2 A 2 ) + \zf >* A,A t {-3Z<V (2 + u 2 A 2 + u 2 A 2 ) 
u e A £ (l - 2u e A e ) + u t A t (l - 2u t A t )j } 



+ A$ 



+ ^^A £ A t 



^(Z e u e A 2 + Z t u t A 2 )-^- 



■ (4 - u e A £ - u t A t + 2u\A\ + 2u 2 t A 2 t) 



P -A,A t 



+ 8 ij 
+ 2 A t 



(Z«) S» (4 + ujA 2 + u 2 A 2 ) + Z™Z™& (ujA 2 + u 2 A 2 ) 



6A £ Z^ + zf^At + \zf > $p 2 AA {zfzf 1 ) $p 2 A t - 4Z<* x >) 



7('j) 







+ -Zf } $p 2 AW ^Zf 'zr'Sp'A 



1 



*»>)]}• 



(3.36) 
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Consider the limiting case in which the derivative terms of the Z' s can be neglected, the above 
expressions can be reduced to: 



(3.37) 



2rf,2 



and 



Zki = Zi + 



A 2 $ 



A 2 $ 2 

1927T 2 



Z t uiAj + N g 1 Z t u t A\ 



(3Z e k 2 + Uj)u t | jV - 1 (3Z.A; 2 + u t )u t 
Z e (Z e k 2 + u £ )* 9 ^(Z^ + u^ 



Zk,t — z t ., 

in agreement with that obtained in [7] for the one-component case. 



(3.38) 



(3.39) 



IV. RENORMALIZATION GROUP FLOW EQUATIONS 

Equation (3.34) gives the contribution of the modes between k < p < A to the blocked 
action in the one-loop independent mode approximation since the systematic feedbacks from the 
high modes to the low ones are neglected. In order to improve upon such approximation, we first 
consider the case when the cutoff is changed infinitesimally from k — > k — Ak, leading to: 



kd k U k 



k 4 



16tt 



In 



1 + - 



A$ 2 /2 



Z t k 2 + 



+ (N - l)ln 



1 + - 



A$ 2 /6 



Z t k 2 + /J 2 



(4.1) 



which is a linear partial differential equation. The equation is not yet suitable for the systematical 
repetition of the elimination of the modes since the right-hand side of (4.1) is derived by using 
the specific potential, (3.21). Since elimination of modes changes the specific structure of the 
lagrangian, it is better to start the whole computation with a general potential. Upon replacing 
the ^-dependent terms on the right hand side of (4.1) by 



u { k 2) m 



d 2 U k m 
dd> 2 ' 



(4.2) 



one obtains a new RG equation: 



kd k u k m = 



k 4 



16tt 2 



In 



ZkAQW + u^Ho) 



+ (N- l)ln 



Zk, t ^)k 2 + Ui 22 \^) 

z k , t (o)k 2 + ul 22 \o) 



(4.3) 
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which accumulates the effects of the eliminated modes in a systematic way as we lower the cutoff. 
Contrary to (4.1), eq. (4.3) is now a nonlinear partial differential flow equation. In the same 
manner, we can write down the corresponding RG flow equations for the Z's: 



kdkZ k ,t = - 



k 4 



16tt 2 



Ah 



- (i#) V(l + tf^Ay] + (N - l)U^Al t {zQur ] (1 " 2U™A k , t ) 



- ^(Jg)A M A M (4 + (^ n) ) 2 A 2 M + (^ 22) ) 2 Aj >t ) + i$A M 



and 



fc 4 f 



*^M<«) = " li^j-^M* 2 [(^ U) ) 2 A 4 M + Aj >t (3 + ^(^ 22) ) 2 A 2 M ) 

-(^) 2 ^fc 2 At(2+(C/f ) ) 2 A^ f ) 

^ -3^> 2 (2 + (^ 11 )) 2 Al > , + (^ 22 )) 2 Ay 

V^AmP - »/f%) + ^ 22) A M (1 " 2^ 22) A M )] } 



(4.4) 



+ -Z$ } $A M A M 



(in) 



+ C/f V A,A t [^ 221) A 2 , + Z k M 22) At t ) 



2 

/c 2 



(4 - ^")A M - U^A k , t + 2{U^fAl, + 2(^) 2 A 2 )t ) 



- — At. P A 



+ s ij 



(Z$) ^•(4 + (^ 11 )) 2 A 2 M + (^ 22 )) 2 A 2 M ) 

+ ig;)iiy* 2 ((^ ii >)^+(^ ) ) 2 Aj >t )" 

(11), ~(fcfc) A . ,°kMk)^2 



+ 2A M + ^^^A M (iii fc j^>^ 2 A fc>t - zf> 



where 



and 



Ok 



i-l(fc 2 ) 2 A M A M (4;^) 2 i 



A-j = Z M ($)fc 2 + ^($) 
A"]=Z M ($)fc 2 + ^ 22) ($). 



(4.5) 
(4.6) 

(4.7) 
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So far the only approximation we made was to make truncation in the derivative expansion 
up to 0(d 2 ). The complicated flow equations can further be simplified if one also truncates in 
the amplitude of the fluctuations which is another expansion parameter of the Landau-Ginsburg 
method. By taking a constant background $o along the longitudinal direction and replacing $o 
by the general inhomogeneous $, we have 



V. 



(22) 



k 

(111) 



k 



221) 



and 



Z 



(i) 



2(U k + 2U k ^ 2 ) 
4{W k +2U' k "<S> 2 )§ 



(4.8) 



zf = 
Zi n) = 4(Z k ^ 2 + Z k ) 
, = 2Z k 8v, 



(4.9) 



where $ points in the 1 direction and the prime denotes differentiation with respect to <E> 2 . The 
three coupled non-linear differential RG flow equations now become (see Appendix B) 



kd k U k ($) 



16tt 2 



In 



+ (N — l)ln 



Z M ($)fc 2 + 2(^($)+2^(^)$ 2 ) 
Z k A0W+2U' k (0) 

Z k , t ($)k 2 +2U k (<S>) 
Z k , t (0)k 2 + 2U k (0) 



(4.10) 



kd k Z ke = 



k l 



16tt 

fc 4 



16tt 



2{ a k,i + a'k,e + a k,£ + a>k,t ~ 2& m) 

4A 2 ^ 2 (-(Z' k /k 2 \l + 4A 2 M (U' k + 2U' k ^) 2 



+ 8Z k , e A% e {U' k + 2U' k ^ 2 ) {3U' k ' + 2U' k "<$> 2 ) 2 

- 2Z k/ (SU' k ' + 2U' k "<f> 2 ) [l - 2{U' k + 2U k ^ 2 )A k , e + 8(U' k + 2U' k ^ 2 ) 2 A 2 k/ 
+ 4(N - 1) A 2 M d> 2 Uu' k A k>t [A M (2Z k , t (U'>) 2 - (Z' k;t ) 2 k 2 U' k ) + Z' Kt U' k \l - 4U' k A k , t )] 



(4.11) 



Z'kAZktk 2 + 2C#) + \2{Z'l^ 2 + Z' kA )A kA + 2{N- l)Z' k/ A k , t , 
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and 



k 4 ( \ 
kd k Z k)t = -y^( a M + a 'k,t + ak,t + a* kjt - 2b k ,t) 

= ~^^A k/ A k ^ 2 i^-Z'^k 2 [z' k>t + 2(U k ) 2 Al t (Z' K£ + 2Z' kjt ) 
+ 16Z k M(U' k r&h - (Zkt) 2 k 2 (l + 2(U' k ) 2 Al t ) + 4Z' k/ A k , t U> k U' k '{l - 4U' k A k , t ) 



- AZ' Kt U' k ' [l - U' k A k , t + 4{U' k ) 2 A 2 ^ + 4U' k ' 

- Z' ktt (l - (U' k + 2U^ 2 )A k/ + 4(U' k + 2U' k ^ 2 ) 2 A 2 k/ ) 
+ 4Z' k/ U' k 'A k/ (U' k + 2U' k ^ 2 ) [l - 4(U' k + 2U'>$ 2 )A k/ 

Z' Kt + 2(Z' k/ + 2Z' ktt ) {U' k + 2U'l&) 2 A\ t 



4Z k/ U> k \U' k + 2U> k ^ 2 )A 2 k/ 



Z'kik 2 



+ 4(Z' k \ t $ 2 + Z' ktt ) A M + 2(N + l)Z' Kt A Kt 



In the limiting case where N = 1 and the derivative couplings are neglected, we have 
kd k Z kA = -^Z fc ^ 2 A 4 M (^ + 2 ^ 2 )( 3 ^ + 2 ^^ 2 ) 2 - 



(4.12) 
(4.13) 



V. ASYMPTOTIC SCALINGS 

a. Ultraviolet regime 

As can be seen from (4.8), the O(N) scalar model possesses a natural mass scale: 

, d 2 



U k=0 {a) = 4<j 2 U' k \ 



(5.1) 



where a =< $ >, the vacuum expectation value. This mass scale is what separates the UV from 
the IR regimes. For sufficiently deeply in the UV or the IR regime, and in the linearizable vicinity 
of the fixed point(s) one finds asymptotic scaling. The calculation presented above reproduces 
the expected perturbative results in the UV scaling regime. To make contact with the usual RG 
equation where only the relevant and the marginal operators associated with the UV fixed point 
are followed, we turn to the RG coefficient functions for the longitudinal modes: 



where 



(3 n (k) = kd k g n , 



g n = dl 1 U k {a). 



(5.2) 
(5.3) 
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In order to find the true critical exponents we have to make the coupling constants dimensionless 
by the help of k. To this end we introduce 



(3 n {k) = kd k g n , 



(5.4) 



with 



~g n = k n - A dl 1 U k (a). 



(5.5) 



The naive power counting for determining the sign of the critical exponents is especially trans- 
parent in this case. In fact, (5.2) now becomes 



p n (k) = k n - 4 (kd k +n-4:)g n , 



(5.6) 



where kd k g n is treated as small in perturbation expansion. Thus, g n is classified as being relevant 
or irrelevant for n < 4 or n > 4, respectively. 

In the case of spontaneous symmetry breaking for which a ^ 0, the evolution equations of 
the quadratic and quartic coupling constants take on the forms 

[3 2 = lm = kd k (4a 2 U k / (a)) 



/c 4 



A M (3Ctf + 12a 2 tC + Aa 2 Ur) ~ 3a 2 A 2 M (3t^ + la^'f 



4tt 2 



+ (N — 1) [A M (Ctf + 2a 2 U> k ") - 4a 2 A 2 ^) 2 ] i, 



(5.7) 



and 



(3 A = kd k [4(3E# + 12a 2 U^' + Aa A U' k "') 
3k 4 



2^ 



jsA^' - 2A 2 fc) ,( 9 (^) 2 + lWotU'W + 224a\U> k "f) 



+ 96A^a 2 (U' k ' + Aa 2 U' k ") (SU' k ' + 2a 2 U'") 2 + 64A 4 ja\3U' k ' + 2a 2 U' k "Y 
+ (N — 1) [a m C#' - 2A 2 M {{U' k r + 12a 2 U' k V' k " + 4a 4 (t/f) 2 ) 

+ 32A 3 M a 2 (^) 2 (^ + 2^') ~ 64A 4 M a 4 (^) 4 l + ■ 



(5.8) 



where the ellipsis denotes 0(U^) terms, and the inverse propagators now become: 

A M = ^k,ek 2 + m 2 R 
Afc j = Z Kt k 2 . 



(5.9) 
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In addition, there are also flow equations for terms which are odd in a due to spontaneous 
breaking of symmetry. For example, we have 

[3 3 = kd k Ua(3U k l + 2a 2 U k 1 ') 



ok 



= ~l 15A M L/f - 18A^(3Ctf + 2a 2 U' k ")(U' k ' + ^ 2 U' k ") + 16A 3 M a 2 W + 2e 2 U' k ") 2 



+ (N - 1)A M 3Ctf' - 6A M Ctf (Ctf + 2a 2 Ctf') + 16A^ 2 (C/£) 



,2 „2( T t"\3 



+ 



(5.10) 

The RG coefficient functions for the wave function renormalization constants are given as: 



= j 4 A 2 ^ 2 {-{Z' k /k 2 [l + mu'<fAU*Y 

- 2Z' kti (3U}; + 2U' k "<j 2 ) \l - 4UgA kfi o* + 32(C/£) 2 A 2 > 



(5.11) 



+ 16Z M ^A 2 ^(3^ + 2^Vf | -4(iV-l)^(^^ 2 +2^)A 2 )t a 2 



+ 12(^> 2 + i^)A M + 2(N- l)i^A M , 



and 



7fc,t = 2 k,t kd k Z kjt 

Z~ x k 4 ( ( 

+ 8^{4Z fc ,,(^) 2 A 2 fc ,,a 2 - Z' ktt [l - UZA^o* + 8(^) 2 A 2 > 2 ) 2 ] } 
+ %Z' k ^U'lf A M <7 2 [l - 8C#A M <r 2 - + 2i£ >t )fc 2 A M <r s 

+ 4(^> 2 + ii >t )A M + 2 (AT + l)l£ >t A M 1. 



(5.12) 



Using = Xr/12 by neglecting 0(U k '), i.e., dropping the irrelevant coupling constants in 
the UV scaling regime, and setting m 2 R = A/jcr 2 /3, the above expressions become: 



02 = 



X R k 4 
16tt 2 



!< 



JV-1. 
Aw H o — A fc,t 



)-A fl a 2 ( 



A 2 I iV - 1 A 2 



))■ 



(5.13) 
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03 



a\ 2 R k 4 
48tt 2 



jsA 2 ^ - 2X R a 2 A kA ) + (N - 1)A 2 M (l - -^a 2 A M ) j, (5.14) 



+ 



(~ - 1) s 



27 



1 - 4A fl a 2 A M + 2\\a 4 A\, 
1 - 12A fl c7 2 A M + 2A|a 4 A 2 5< 



(5.15) 



7M 



^k,£ kdkZk, 



A 2 . 2 , 2 ,21 A 3 



+ — Z k,iA 2 kl o 2 



12 



ki 



A 



.2 , 2A 2 2 •_, 2 



1 - -A M a 2 + ^-KM 



4(N-l)Z' k /Z' ktt k 2 + -)Al t a 



(5.16) 



+ 12(i£> 2 + Z' k , e )A k , e + 2{N- l)Z' k/ A k , t 



7k,t = Z k I kd k Z, 



k,t 



^k,t 

16tt 2 



2A kj gA kyt a 2 j —Z' kt k 2 (2Z' kl + Z' k t j 



2A 



A 



A 5 



-Z M A 2 ^ 2 - Z' ktt {l - -A kJ a 2 + -Al/a 2 ) ) 
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A 2 r 2A 

+ j^Zk,e A k,t<T 2 [l - y A M a 2 - + 2Z^ t )k 2 A k ^ 



(5.17) 



+ 4(^> 2 + Z' Kt )A k/ + 2(N + l)Z kt A k j 



The leading order part of the above coefficient functions for k 2 jm 2 R 
l + 0(X R ): 



02 
03 
04 

lk,e = 
lk,t = 



N + 2 



48tt 2 



\ R k 2 



aX 



R 



48. 2 ^ + 8 ) 

Alt (N + 8, 

16n 2[ 27 ' 

A| /^ 2 x2 
487r 2 I k 2 ) 



48tt 2 v 3/c 2 



Vql2/ 



0. 



cx) and Z fe Z, 



k,t 



(5.18) 
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The last two lines are in agreement with the well known fact that there is no need for wave 
function renormalization at the one-loop order. In fact, their coefficient functions are vanishing 
in the ultraviolet. When symmetry breaking appears, Zk,i and Z^t evolve differently from each 
other as we leave the UV regime. 

If one wishes to compare the above results with the usual e expansion in critical \(p 4 theory 
in 4 — e dimension, it is sufficient to consider a symmetrical theory in which one sets U k (0) = 
and defines a dimensionless coupling constant = 12k~ e U%(0) = 12k~ e \k- From (5.8), one is 
led to 

. „ r r 3A 2 fc r 1 JV-1 1 i 3U"'(0)k 2 i 5 iV-li 

which for Z^i = Zk,t = 1 reproduces the leading order results obtained by standard perturbation 
method. 



b. Infrared regime 

The scaling is more involved in the IR regime where the expansion k 2 « m 2 R should be 
applied [6]. The strong non-linearities of the RG equation prevented us from constructing the 
corresponding scaling operators. Instead we only argue that relevant operators must exist in the 
IR scaling regime. 

Let us first begin with the naive argument outlined in the Introduction where Z is set to be 
unity. In fact, due to Goldstone's theorem which asserts U^Q(a) = U' k=0 (a) = 0, we take 

Ul 22 \a) = 0(k?) (5.20) 

with p = 2. In this limit with U^ 2 \a) = ck 2 , the leading order contribution to the (3 function 
becomes 



1} i Uk {a) ^ (l + 0(fc 2 )), 



(5.21) 



16tt 2 V. (l + c)fc 



where 



rr (i22)/ , d U k {cr) . TT „ Xr /r ^ 

u * (a) = a$p$7 = 4at/fc = T (J + ---- (5 ' 22) 

( 122) 

Since Ul (cr) is finite for non-critical system inside the symmetry broken phase the j3 functions 

( 122) 

develop power-like IR divergences when W k '(a) > 0. 

This simple argument relies on the leading order of the derivative expansion. Having gone 
through the computation of the wave function renormalization constants we can verify that our 
conclusion remains valid in the next order of the derivative expansion, too. In particular, we 
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show that the IR power singularities of the (3 functions persist for the RG equation (4. 10) (4. 12). 
In the presence of a nontrivial wave function renormalization constant the (3 function reads as 



(-l) n k n fZ { ^Aa)k 2 + ui 122 \a)\ n 
H) " V Z kit {a)k 2 + U k » J 



(-1)- ( uj 122 \<j) 

16tt 2 1 k x ~ v 



(l + 0(k 2 )) 



where the asymptotic behavior 

Z k , t = 0[k~ v ) (5.24) 

v > was assumed. 

In order to show that v < 1, which would imply the persistence of IR singularity in (5.23), 
we verify the consistency of (4.11) and (4.12) by assuming (5.20), (5.24) and 

Z k ,t = 0(k-»). (5.25) 

A similar power-law dependence on k also applies to the Z' k s and the Z' k n s since differentiation 
with respect to $ 2 does not affect the power counting in k. While Ak,e remains finite as k — > 
due to the presence of mass gap, the transverse propagator is scaled as 

A M = 0(k- 2 ). (5.26) 

Substituting the above scalings in (4.11) and (4.12) for matching the leading singularities gives 

-p = minj 6 - 2p, 4 - p, -4 + 3z/ + p, 2 - p + u, -2 + 2v + p, ■ ■ •}, (5.27) 

and 

— v = min< 4 — p, 4 — z/, — 2^ + 3z/ + 2p, — 2 + 2z/ + p, — p + 2z/ + p, 2, 

1 (5.28) 
-2-p + 3z/ + 2p, v + p, 2 + p, 2-p + z/, 4-2p + z/H L 

Relying on the method of independent-mode approximation in [6] , we have p = 2 which im- 
plies [i = 2 and z/ = 0. Thus, we conclude that presence of IR singularity in the (3 function (5.23) 
remains unchanged after taking into consideration the wavefunction renormalization constants. 
Note that the longitudinal wavefunction renormalization constant Zn is found to be quadratically 
divergent in the IR limit. This is in contradiction with the usual assumption [9] 

< Z < 1, (5.29) 

made in Minkowski space-time. Thus the invariance of Z under the Wick rotation, the one-loop 
RG equation, i.e. ( ~ and (5.29) are inconsistent. It is not clear how to describe the vacuum 
but a deviation from the weakly coupled perturbative scenario is expected due to the Goldstone 
modes. 
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We emphasize again that the difference between our (3 functions (5.7), (5.8) and the usual 
ones in (5.18) stem from the terms 0(a 2 /k 2 ) which are neglected at the UV fixed point. These 
are just the pieces which make the UV and the IR scaling laws different. 

In the course of investigating the IR scaling behavior, our computations are free of IR 
divergences. This is because for a finite value of the adjustable cutoff, k, the elimination of the 
degrees of freedom in (3.6) contains the integration for modes with finite momentum k < p < A 
and the usual singularities of the massless theories do not show up. They appear only when 
k — > 0, the limit where singularity appears in the (3 functions. The determination of IR scaling 
operators is rather involved since it requires a complete resummation of the singularities which 
emerge in that limit. 



VI. SUMMARY 

Most of the applications of the RG are related to the investigation of the impact of local field 
operators which are introduced into the theory in the UV regime. The concept of universality 
which is being supported by the linearized RG equations can be used as long as the phenomena we 
are interested lie within the linearizability of the UV fixed point. Such phenomena can therefore 
be described by a simple hamiltonian containing few relevant operators. This certainly is the 
case for the critical, i.e. massless <^> 4 model near four dimensions where both the UV and the IR 
fixed points are Gaussian. The classification of the scaling operators is valid for all length scales 
in this model. 

The situation is radically different for massive theories. There we have two energy scales, 
the UV cutoff and the mass gap. There is no reason to expect the same scaling laws at both 
sides of the mass gap. In fact, the renormalized trajectory runs towards large values of the mass 
squared as we approach the IR regime. As the vacuum expectation value of the field reaches 
0(l/e), where e = 4 — d then the three-point vertex becomes of order 1 in the symmetry broken 
phase and the expansion around the Gaussian fixed point is not applicable. 

This is the generic situation for high-energy physics where the renormalized trajectory passes 
by the vicinities of different UV fixed points and we find different scaling laws at different energy 
ranges. The usual concept of universality is not applicable here since the relevant coupling 
constants of an UV fixed point parametrize only the physics of a given energy range. The 
trajectory is driven away from the region of the linearizability by the relevant or the irrelevant 
coupling constants as we move towards lower or higher energies, respectively. For the sake of 
definiteness we considered the O(N) model in this paper which has a single finite scale and 
exhibits only two fixed points, an UV and an IR. 

It is usually claimed that there is only one (completely trivial) relevant coupling constant at 
the IR fixed point, the mass. But this claim ignores the IR divergences of the massless theories 
which may generate relevant operators as the observational scale approaches the IR regime. 
Another class of models where the IR scaling might be rather nontrivial is where a symmetry is 
broken spontaneously. This case is interesting because it emphasizes the importance of 

2"7T K 

C= Lk = V (fU) 

the ratio of the observational length scale, k = 2nk~ 1 , and the IR cutoff, L. For L close to 
the characteristic mass scale of the theory there is no symmetry breaking and the evolution 
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equation for the effective coupling constants reflects the symmetrical dynamics. For large but 
finite L the symmetry is still preserved but the spontaneous symmetry breaking scheme becomes a 
good approximation for the dynamics. Spontaneous breakdown of symmetry can take place only 
asymptotically in the limit L — > oo. Since k can never exceed L for a finite system, we have £ < 1. 
Therefore, pattern of symmetry breaking can be uncovered in the evolution equation only for 
£ ~ 0. In fact, for £ ~ 1, one detects the symmetry- restoring long-range slow fluctuations, and the 
observables at this energy scale truly reflect the symmetrical dynamics. Thus the characteristic 
size of the system ~ L should be much larger than the observational length scale k in order to 
recover the usual picture of symmetry breaking. 

There is another rather technical reason for staying in the region of small £. The higher loop 
contributions to the RG equations are suppressed by the inverse of the number of the modes in the 
blocked system which is 0(£ d ). Thus the studies of systems undergoing spontaneous symmetry 
breaking using the one-loop RG equation requires the removal of the IR cutoff by sending L — > oo 
before k — ■> oo, thereby making £ = 0. 

The locality is lost at the IR fixed point which may lead to the breakdown of the derivative 
expansion and global scaling operators. All we know is that classical physics is recovered at the 
IR fixed point of massless theories. In this case the soft particle emission allows the spread of 
the energy from the microscopical to macroscopical length scales. On the other hand, when the 
theory possesses a mass gap, the energy cannot be distributed to arbitrary long distances and 
the IR physics is still controlled by coherent quantum effects, e.g., superconductivity. 

We found IR singularities in the one-loop (3 function for the odd vertices of the 0(N) model 
when the first two orders in the derivative expansion of the renormalized lagrangian were retained. 
This supports the notion of strong coupling IR physics of the Goldstone modes. The amplification 
of the effective coupling strengths can be understood by recalling that the "restoring force" for 
the fluctuations, i.e., the eigenvalue of the small fluctuation operator is vanishing in the IR limit 
of a massless theory. Consequently large fluctuations are always present in the IR regime and 
invalidate the expansion methods. 

The limitation of the concept of universality due to the existence of several fixed points can 
be nicely demonstrated in the O(N) model. Consider the coupling constant of an odd vertex 
in the spontaneously broken phase. Being irrelevant in the vicinity of the UV fixed point, it 
decreases as we move in the IR direction. Its value should be small when we reach the crossover 
region between the UV and IR scaling, at the mass gap. Universality, i.e., the insensitivity on 
the irrelevant initial conditions of the renormalized trajectory seems to be holding down to this 
energy scale. But as we continue our journey in the IR direction our coupling constant starts to 
grow. Although we could ignore it on the UV side of the mass gap, it plays an important role on 
the IR end. Furthermore, its actual value may depend strongly on the ultraviolet initial value 
of the renormalized trajectory. The suppression which produced the universal behavior down 
to the mass gap turns out to be an amplification in the IR scaling regime and the UV value of 
this coupling constant influences the long-distance features of the model. This possibility raises 
questions on the sufficiency of renormalized field theories in describing low energy phenomena. 

We believe that only few nonrenormalizable operators become important in this manner. To 
demonstrate this point consider the Standard Model. Despite all complications in the IR regime 
the experiments performed in the vicinity of the crossover, O(100GeU), can be parametrized 
by the help of the renormalizable coupling constants. Where then is the room for the possible 
violation of universality ? The conjectured strong coupling physics in the IR plays an important 
part in forming the vacuum but remains virtually invisible at higher energy. The only important 
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parameters they provide are the values of the condensates. These condensates appear under the 
disguise of renormalizable coupling constants, say lepton masses in the usual scheme. But the 
relation between the mass and the condensate bridges the energy scale of the mass and zero and 
thus its tree-level form is highly questionable. On the one hand, the order parameters of the 
spontaneously broken symmetries are formed in the asymptotical IR regime. On the other hand, 
they parametrize the effective vertices at the crossover. 

It is reasonable to expect that the only impact of the IR modes on the physics of the higher 
energy processes is the generation of the symmetry breaking condensates. Thus universality 
actually holds when the physics is parametrized by the help of the renormalizable coupling 
constants and the condensates. It is useless if we make an attempt to derive the values of the 
condensates starting with the UV parameters, from one fixed point only. 

This scenario leaves universality unharmed for ferromagnets. In fact, the condensate of the 
nonlinear a- model is a unit vector and there is no possibility of changing its length. In contrast of 
this situation, the physics of the superconductors may show non-universal features. In particular, 
the supercurrent density might depend on non-renormalizable coupling constants of QED which 
are provided by theories of higher energy scale, such as the Standard Model. In turn, the Higgs 
condensate of the Standard Model is a non-universal function of the bare parameters of a GUT, 
etc. 

One would object the speculations about relevant operators for the IR fixed point of a 
superconductor since there is no gap in the physical spectrum. But the massless excitations are 
present in the gauge-dependent sector where the Higgs mechanism relegates them and continues 
to influence the dynamics of the gauge invariant modes. Their presence can be seen from the 
long range confining forces acting between two magnetic charges. In fact, in the absence of 
massless modes all interactions are screened. The situation is similar to the QCD vacuum where 
the long-range confining modes coexist with finite range the Yukawa forces due to the massive 
glueball exchanges. 

In closing we repeat again that our results rely on the derivative expansion. It would be 
of key importance to support or disclaim its validity for the four-dimensional models in the IR 
regime. 
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APPENDIX A 

We collect here the commutation relations which were used in deducing the wavefunction 
renormalization constants. All the relations derived in this Appendix are based upon the simple 
rule 

[ 0, Vll ) = -id^. (A.l) 

Note that the appearance of the negative sign which differs from the conventional definition 
is due to the fact that the space-time traces are evaluated in the plane-wave basis with all p 
dependences being moved to the left of ^-dependent field operators. Repeated use of the above 
gives 

[ 0, p 2 ] = -d 2 - 2ip^4> (A.2) 
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i, P»P 2 ] = - {P»d 2 4> + 2p v d„dp$) - i {p 2 d^ + 2 PflPl/ d v $) + ■■■ (A. 3) 

[ 4>, (P 2 f ] = -3p 2 d 2 4> - Aip 2 p^ +■■■ (A.4) 
[[4>, p 2 ], p 2 ] = -Ap^d^ +■■■ (A.5) 



(A.6) 



>, A] = ZA 2 [p 2 , <p]+Z 2 A 3 [p 2 , [p 2 , 0]] + -.. 
= ZwA 3 <9 2 + 2iZA 2 p fl d fl 4> +■■■ 

[ O b , ] = [ A ] = (A.7) 



[ 4>, p 2 A ] = -Z- X u\ <j>, A } = -u 2 A 3 d 2 4> - 2iuA 2 p fl d fl <fi +■■■ (A.8) 

[ 4>, A Pfl ] = ZA 2 (uA PtI d 2 4> + 2p v d v dJ>) + i A (2ZA PflPv d v 4> - d^) + ■■■ ( A.9) 

[ 4>, (p 2 fA ]=Z-'[^P 2 ]+ Z~ 2 u 2 [ 0, A ] 

= -p 2 A{l + uA + u 2 A 2 ) d 2 cj)-2ip 2 A(l + u A)p^d^+ ■■■ 

[ 4>, p 2 A Pfl } = -uA 2 {uAp^d 2 j) + 2p v d v d ll $) - iA(2uAp lx p v d v ^ + p 2 d fl 4>) +••■. (A.ll) 

[ 4>, A e A t ] = A e A t {(Z t u t A 2 t + Z e u e A 2 - Z £ Z t A e A t p 2 )d 2 4> + 2i(Z e A e + Z t A t )p^d^} + ■■■ 

(A.12) 

[4>, A 2 ] = —ZA 3 (l — 3wA) d 2 (f) + 4iZA 3 p ll d fl (f> + • • • (A.13) 

[ 0, p 2 A^A t ] = A £ A t \(Z t p 2 u t A 2 + Z £ p 2 uiA 2 e - Ul u t A £ A t )d 2 ^ 

, ^ ~1 (A ' 14) 
+ 2i(l - it*A* - u t A t )p^d^(j) \ H 

[0, p 2 A 2 ] = «A 3 (2-3wA)d 2 + 2zA 2 (l-2wA)p M d M 0+--- (A.15) 

After generating the derivative terms with the above commutation relations, one may un- 
tangle the x— and p-dependent terms with the following useful relations: 

(p 2 h - 2ip^h - d 2 f 2 + f 3 )A = A[p 2 h - d 2 f 2 + f 3 + ZA[p 2 (l + uA)d 2 h + uAd 2 f 3 ] 

+ 2ip fl A[-ud^f 1 + Zd fl f s ]} + ---, 

(A. 16) 
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(A.17) 



(A.18) 



{p 2 h - 2ip fl dJ l - d 2 f 2 + h)A(p 2 gi - 2ip v d„ 9l - d 2 g 2 + g 3 ) 
= AL 2 (p 2 h + f 3 ) 9l + {p 2 h + f 3 )g 3 - {p 2 h + f 3 )d 2 g 2 - (p 2 9l + g 3 )d 2 f 2 
+ [p 2 (2 - u 2 A 2 ) 9l + Zp 2 A(l + uA)g 3 ] d 2 h + [(1 + uA - u 2 A 2 ) 9l + ZuA 2 g 3 ] d 2 f 3 
+ 2tp II A(p 2 g 1 + g 3 ) {-ud^h + ZdJ 3 ) + ■ 
which for g = f reduces to 

{p 2 h - 2ip^dJ 1 - d 2 f 2 + f 3 )A{p 2 h - 2ip v d v f 1 - d 2 f 2 + / 3 ) 

= a|p 2 (p 2 /i + /s)/i + {P 2 fi + fa) fa ~ 2(p 2 /i + h)d 2 f 2 
+ p 2 {2- u 2 A 2 )hd 2 h + ZA 2 uf 3 d 2 f 3 + (2 + uA- 2u 2 A 2 )f 3 d 2 f 1 
+ 2ip^A{p 2 h + / 3 ) {—ud^fi + Zd^h) + 

Similarly, we have: 

{p 2 h - 2ip^fr - d 2 f 2 + / 3 ) (p 2 9l - 2ip v d v9l - d 2 g 2 + g 3 ) = (p 2 ff l9l + p 2 (f l93 + f 39l ) 
+ f 393 - 2ip^p 2 9l + <7 3 )<Vi + (p 2 h + / 3 ) {d 2 9l - d 2 g 2 ) - (p 2 9l + g 3 )d 2 f 2 + • • • , 

(A.19) 

and 

(fi - d 2 f 2 - 2ip li d li f 3 ) [g + (p 2 fgi + p 2 g 2 - d 2 g 3 - 2ip 2 p„d u g A - 2ip u d v g b ] 

= fm + (P 2 ) 2 fi9i +P 2 fi92 - d 2 f 2 g -p 2 {3d 2 f 1+ p 2 d 2 f 2 + 2p 2 d 2 f 3 ) 

- (d 2 f 1+ p 2 d 2 f 2 +p 2 d 2 f 3 )g 2 - hd 2 g 3 + p 2 (3f 1+p 2 f 3 )d 2 g 4 + {2f x +p 2 f 3 )d 2 g 5 + ■■■. 

(A.20) 



APPENDIX B 

In computing the effective blocked action Sk, one encounters the N x N matrix M of the 



form: 



where 



and 



a T =(a e zf } ... a lZ \ N) ) a^\p 2 ^A t 



b = 



( a t zP \ 



\a t zp; 



a t = ^p 2 $ A t . 



(B.l) 
(B.2) 

(B.3) 



27 



One can easily verify that its inverse M 1 takes on the form: 

6 -6a T 



M~ = 



where 



Employing the relations: 



-9b I + 6ba T ' 
v/a- 1 _ ri ~„ /^hV 1 



9={l-a T b) =[l-a £ a t (Zn ] 



oo , ^ \n+l 



n=l 



and 

oo 

^(-l) n Tr' 

n=0 

we readily obtain 



)n 
KqHk 



= Tr' 



+K- i dK Q y 



K~HK, 



UJ1 a 



(B.4) 
(B.5) 

(B.6) 
(B.7) 



5S 



I = -Tr'ln^o + SK + 8K X + SK 2 ^ = -Tr'ln^o + 6K } - -Tr' {k^^K^-^K^ 

- oo oo 

1 r / \ n -i 1 r / \ n 



n=0 



n=0 



= ^Tr'ln^o + ^o) - [Kq-^SK^Kq^SK^ + ^Tr' [(l + K^SK^j 1 K^ l 5K x 
+ ^Tr' [(l + K^SKo) "Vo" 1 ^] + • • • , 



In terms of matrix elements, 

Tr'[(l + K^SKoy^^SKa] = 



(B.8) 




#A,(5K a ) n + 6A t a t a t zyzy(5K a y i 



x Jp 



(B.9) 



+ A t - 0(a*A t + a t A^)zf (<5K Q ) a 



and 



(B.10) 



+ 



A^K^At^y + AtiSK^AeiSK^ 



Nil 



a? 



The complicated commutator algebra can be simplified by noting that the matrix elements take 
on the forms: 



(5-ftTi) ab = p 2 ax - 2ip fl d fI a 1 - d 2 a 2 + a 3 
(SK 2 ) ab = p 2 fofr - fod 2 4> c + fa, 

28 



(B.11) 
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ai = ~zf ) ft + zf c 



a 2 = \z<pp 
as = A^o^ 1 



+ Z, 



a 1 = \{zf^ + ~Z ( t 1) ~^ + Z^ ) ^) 



«2 = '<j>~ + Z\ 



zf- 1 



1 

2" 



+ ^zf ^ 



(3, = z^p + \z\ cd) ^ d 



{SK^f : I 



(SK 2 ) 



li 



02 



A 



A ^(zf^ + zf^ 1 ) 



' (3[ = -{z[ ic) ft + Z? c) fi + Z ( t cd) d^4> d ) 
2 



(SK 2 ) 



Substituting (B.14) - (B.16) into (B.9) for 5K 2 yields: 

Tr'^l + Ko-^Ko)" 1 ^" 1 ^] = J J {b 11 4> 1 d 2 4> 1 +b lJ 4> t d 2 ft} + - 



where 



6 n = -i(3zi n) (9A, + zf } A f ) - |z« <& A,A fP 2 [zf zf ^ A t p 2 - 8Zf >] 

o 



6« = "o^ 



0A,Z t (11) + Zf fc) A t + -Zf ^ p 2 A,A,(zf Zf >*Qp 2 A t - 4zf 1} ) 



zf > + ^zf $ oP 2 A^zf >Z^>$ P 2 A t - Zp' 1 )) 



4 
1 



(B.12) 



(B.13) 



(B.14) 



(B.15) 



(B.16) 



(B.17) 



(B.18) 



29 



For the 5K\- dependent terms, after much tedious algebra with the help of the relations 
found in Appendix A, we have: 

AeiSK^AeiSKi) 11 = A 2 £ ^{p 2 fa\ + 2p 2 a x a 3 + a\ + 2p 2 a 1 (d 2 a 1 - d 2 a 2 ) - 2a 3 d 2 a 2 
- p 2 u 2 Aja 1 d 2 a 1 + Zu £ A 2 £ a 3 d 2 a 3 + (2 + u £ A £ - 2u 2 A 2 £ )a 3 d 2 a 1 H 1 



a n = -A 2 UZ^) p 2 (l + u 2 A 2 ) + Z\ L> \<$> {1 - u £ A £ + 2u 2 A 2 ) - Z £ u £ A 2 X 2 $ 2 

aij = -Zfzfy 2 u\A\, 

AtiSK^AtiSK^ = a' n <^ W + a'^d 2 ^ + • • • , 

a' n = (N- l)«tA?{A t [^*gZ t - (Z t (1) )V«t] + ^o^ (1) (l - 2u t A t )} 



-zf^ A?(^ ( V + ^o) 



(B.19) 
(B.20) 
(B.21) 



(B.22) 



(B.23) 



~a 11 = -](Z^) 2 p 2 A £ A t (2 + u 2 A 2 ) 

a ij =A £ A t { \zf j) ^ Q \-^z{ 1) p 2 {l + u 2 t A 2 t ) + \^ Q u t A t {l-2u t A t ) 



-5 lJ 



6 

\{Z?>) V(2 + u 2 A 2 ) - (^oW? - ^(2 - u t A t + 2u 2 A 2 \ 



]z™Z^ 2 p 2 u 2 A 2 }, 



A t (SK 1 ) il A £ (SK 1 ) li = a^d^+a^d 2 ^, 



(B.24) 



(B.25) 



a 



l 1 = -l(zP) 2 p 2 A £ A t (2 + u 2 A 2 ) 



1 



A £ A t < -zy j ><$> 



-3Z t (1) p 2 (1 + u^) + A$ «<A*(1 - 2u £ A £ ) 

(i) 



\{Z { t 1) )p 2 {2 + u 2 A 2 ) - (^o^A 2 - ^-(2 - u,A, + 2«jAj)) 



(B.26) 
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In arriving at the above expressions, 0(4) invariance has been used: 



/ 0(9^0 / p^p v = - I (j)d 2 (j) I p 2 . 

J x J v J x J v 



(B.27) 



The coefficients obtained above become much simpler when O(N) symmetry is invoked which 
allows us to make the following substitutions when seeking for the RG flow equations: 



and 



U£ 

^ 3 



r 

#11) 



ui 11) = 2(U' k + 2U'^ 2 ) 



U { k 111} =4(3U^ + 2U^$ 2 )$ 



(B.28) 



= 2Z k $ 
zf = 



(B.29) 



where the prime notation denotes differentiation with respect to <£> 2 and we have chosen the 
transverse component to be along the 2-direction. Thus we have: 

6 M = -6{Z' k \ e $ 2 + Z' k ^)A k/ - (N - l)Z' ktl A ktt 

6 M = -2{Z'l^ 2 + Z' kit )A kii - {N + l)li |t A M 



(B.30) 



a k , e = -4A 2 ^ 2 \ (Z' k /k 2 [l + 4A 2 M (U' k + 2U' k ^ 2 ) 2 ] + 2 (3U' k ' + 2U' k "$ 2 ) 

1 - 2(U' k + 2U' k ^ 2 )A k/ + 8(U' k + 2V' k ^ 2 ) 2 A\, 



x Z 



(B.31) 



4Z k , e Al £ (U k + 2U' k ^ 2 ) {3U' k ' + 2U' k "&) 



o' M = 4(N - l)A 2 ^ 2 ^-Z' k/ {Z' kjt k 2 + 2UH) 

+ 4U k A k , t [A M {2Z Kt (U' k ') 2 - (Z^ t ) 2 k 2 U' k ) + Z' Kt U' k \l - 4C#A M ) 

~a k , t = 2A M A M d> 2 1 16Z k , t U' k (U'k) 2 A 2 M - Z' k/ k 2 [z' Kt + 2{U' k f A 2 M (Z' k/ + 2Z' Kt ) 
- (Z' k ^) 2 k 2 (l + 2(U' k ) 2 A\ t ) + 4Z' k/ A k ,tU' k U' k ' (1 " 4UL A M ) 
-4ZL t U'Al-U' k A Kt + 4{U' k ) 2 Al 



(B.32) 



(B.33) 
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at tt = 2A M A M d> 2 AU' k ' AZ Ki U' k \U' k + 2U' k ^ 2 ) A\ t 



and 



- Z' ktt (1 - (U' k + 2U' k ^ 2 )A k/ + A(U' k + 2U' k ^ 2 fAl 
+ AZ' k ^A K , (U' k + 2U' k '&) [l - 4(U' k + 2U' k '$ 2 )A k/ 

Z' ktt + 2(Z' ktt + 2Z' ktt ) {U' k + 2U' k % 2 ) 2 Al 



~ Z'k,ik 2 



a k,t — a' kt — a kj £ — a* k£ — 0. 



(B.34) 



(B.35) 
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